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Abstract 



We have studied the lightest masses in the fermionic sector of an anomalous U{iy extension of the Minimal Supersym- 
. . . jnetric Standard Model (MSSM) inspired by brane constructions. The LSP of this model is an XWIMP (extremely weak 
interaction particle). We have studied its relic density in the cases in which there is mixing with the neutralinos of the 
MSSM and in the case in which there is not such mixing. We have showed that this extended model can satisfy the 
WMAP data. To perform this calculation we have modified the DarkSUSY package. 



1. Introduction 



There has been much work recently to conceive an in- 
3 Itersecting brane model with the gauge and matter content 
rrnof the Standard Model (SM) of particle physics dSSIl. 
'"^ 'One of the features of such models is the presence of extra 
f) anomalous gauge C/(l)'s whose anomaly is cancelled via 
On.the Green-Schwarz mechanism. The anomalies of these 
^ models are cancelled without assumption on the newly in- 
'troduced quantum numbers. One model of this type have 
been introduced and discussed in 0] and 0. We have 
[studied the compatibility of this model with the WMAP 
data (see also 13, HI for related work). In this model 
we have two new contribution to the neutralinos mass ma- 
trix: one coming from the superpartner of the Stiickelberg 
boson (Stiickelino) and the other from the superpartner 
of the gauge boson mediating the extra U{\)' (Primeino). 
By taking some simplifying and reasonable assumptions 
on the fermion masses entering the soft supersymmetry 
breaking lagrangian, the LSP turns out to be mostly a 
•mixture of Stiickelino and Primeino. It is also easy to re- 
alize that, given the simplifying assumptions mentioned 
52^ above, the LSP is interacting with the MSSM particles 
& with a coupling suppressed by the inverse of the mass of 

24 the extra gauge boson of the theory which must be at least 

25 of the order of the TeV for phenomenological reasons. The 

26 LSP is then an XWIMP (Extra Weakly Interacting Mas- 

27 sive Particle), a class of particles which have already been 

28 studied in literature [ij, The cross section of these 

29 LSPs is too weak to give the right relic abundance. This 

30 is why one has to resort to coannihilations with NLSPs. 

31 In our case, the cross section for the annihilations of the 

32 LSP with the NLSP and that for the coannihilations of 

33 the two species differ for some orders of magnitude. Once 



in 
It 



again this situation is not new in literature IJ, Il5l . |16 | 
but needs to be treated carefully: the two species will not 
decouple as far as there will be some MSSM particles to 
keep them in equilibrium. Moreover these particles must 



38 be relativistic so that their abundance is enough to foster 

39 the reaction. There are two cases. First, the anomalous 

40 sector of the neutralinos mass matrix can be decoupled 

41 from the MSSM sector: this implies that the LSP can be 

42 either pure anomalous (Stiickelino-Primeino mix) or pure 

43 MSSM. This is called the no mixing case and we have 

44 studied it in the case in which the LSP is pure anomalous 

45 and coannihilates with a pure MSSM NLSP. The second 

46 case, called the mixing case, is that in which there are mix- 

47 ing terms between the two sectors of the neutralinos mass 

48 matrix. So we can have LSP and NLSP that can be a gen- 

49 eral mixing of the six gauge eigenstates of the theory. We 

50 will show that also in this general case the mixing is con- 

51 strained and then we will study the case in wich the LSP 

52 is mostly anomalous with fractions of MSSM eigenstates, 

53 while the NLSP is mostly of MSSM origin whit small frac- 

54 tions of anomalous eigenstates. 

55 In the subsequent sections we will describe the results for 

56 the relic density in both cases. We have obtained that 

57 there are regions of the parameter space in which the 

58 WMAP data can be satisfied in both cases. These re- 

59 suits are obtained using a modified version of the Dark- 

60 SUSY package which takes into account the couplings of 

61 our model. 

62 2. Model 

63 Our model 'i] is an extension of the MSSM with an 

64 extra U{1). The charges of the matter fields with respect 

65 to the symmetry groups are given in table 1. 

66 Gauge invariance imposes constraints that leave inde- 

67 pendent only three of these charges. In this model we have 

68 chosen: Qq, Ql and Qh^- The anomalies induced by 

69 this extension are cancelled by the GS mechanism. Each 

70 anomalous triangle diagram is parametrized by a coeffi- 

71 cient &2°'' (entering the lagrangian) with the assigment: 
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The mass of the extra boson is parametrized by Afy(o) = 
^ 1 T'eF, where 50 is the couphng of the extra 
[/(!)'. The terms of the Lagrangian that wiU contribute 
to our calculation are |6|, 181] : 



Lstuckolino = ^V'SCT^ 9^-05 " V2b3^s>^'-"^ 

-^E4'^'Tr(A("V'^a^FW)^^ (6) 

h.c. 



^(0 



h('^a''a''Fl°}i:s + (0 O 1) 



The parameters fej ^-re inversely proportional to My(a) , so , 
they are much smaller than the couplings of the SM. This , 
implies that the primeino and the stiickelino are XWIMPS. , 



The neutral mixing matrix is: 



Defining an = go Qh^ 




(7) 
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where Cw is cos{9w), sw is sin{0w)- gi, 52 are the cou- 
phngs of the SM electro-weak SU{2) x U{1) group. The 
structure of this matrix leaves the electromagnetic sector 
and the related quantum numbers unchanged with respect 
to the MSSM ones. We can see from eq. ^ that the case 
of mixing between the anomalous sector and the MSSM 
sector is that in which an ^ 0, while the case of no mixing 



is that in which an = 0. From the definition of an it's 
clear that an — -^Qh^ — 0. 

Last, we remember the general form of the neutrali- 
nos mass matrix at tree level: 
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where Ms and Mg are the soft masses of the stiickelino 
and of the primeino, respectively and Vm Vd are the vevs 
(vacuum expectation values) of the Higgs fields. We can 
see that in case of no mixing {Qh^ = 0), there are no 
terms different from except for the two diagonal boxes 
that refer to the anomalous and the MSSM sectors. 

3. Boltzmann Equation and coannihilations 

The relic density of a certain particle is given by the 
Boltzmann Equation (BE) flTll : 



dn 



an -^—^ 

— = -3Hn- y \UtjVij){n^r 



tr,,- is the annihilation cross section: 



ff 



n^nf) (9) 



(10) 



is the absolute value of relative velocity between i-th 



and j-th particle, defined by: 



is the 



number of particle per unit of volume of the i-th specie 
(index 1 refers to LSP), n^'^ is the number of particle per 
unit of volume of the i-th specie in thermal equilibrium. 
n — ni, H is the Hubble constant and v is the relative 
velocity of the initial particles. 



If the LSP is an XWIMP, it is known [12| that its relic den- 



sity does not satisfy the WMAP observations [18[. So, in 
order to achieve the sperimental results, we need coanni- 
hilations. This phenomenon occurs when there is at least 
one particle with mass of the order of the mass of the 
XWIMP LSP. If the interaction of the coannihilating par- 
ticles are strong enough they can lower the relic density 
of the LSP by several orders of magnitude and so we can 
obtain results in agreement with the sperimental data. 
Although a general tractation can be found in [l^ , for our 
purposes we are interested only in the cases in which the 
number N of the coannihilating particles is 2 or 3. In the 
case = 2 we have the following thermal averaged cross 
section: 
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where Q = /n^. The first term in the numerator can 
be neglected because the Stiickehno annihilation cross sec- 
tion is suppressed by a factor {h^2^Y with respect to the 
MSSM neutralino annihilations (see the previous section) 
and thus {cfhv) <C {(^22^). 

With the same assumptions of the case N ^ 2 the formula 
for the case = 3 is: 



{<T22V)QI + 2(cr23v)Q2Q3 + {(^33^) 



with 



(I + Q2+Q3) 



eq* 

\ 

eq 
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This expressions show that the factors Qi can raise/lower 
the thermal averaged cross section and thus modify the 
related relic density. However to obtain the result for the 
relic density it is necessary to numerically solve the BE. To 
achieve this we have used the DarkSUSY package, adding 
the interactions introduced by our extended model. In the 
following section we describe our results. 



124 4. Numerical computations and results 

125 To perform a numerical analysis on the extended model 

126 we have modified the DarkSUSY package adding new fields 

127 and interactions introduced by the anomalous extension. 

128 The free parameters that we use in our numerical simula- 

129 tions are the seven ones used in the MSSM-7 model: the n 

130 mass, the wino soft mass M2, the parity-odd Higgs mass 

131 Maq, tgP, the sfermion mass scale msf, the two Yukawas 

132 at and a;,. We add to this set five parameters which de- 

133 fine the U{\)' extension: the stiickelino soft mass Ms, the 

134 primeino soft mass Mq, the U{1)' charges Qh^, Qq, Ql- 



4.1. No-mixing case 

As we have seen if there is no-mixing we have to im- 
pose = 0. From eq. this implies that the LSP 
can have either pure anomalous origin or pure MSSM ori- 
gin. We are interested in the first case. Furthermore, we 
assume for simplicity that the LSP is a pure stiickelino. 
Because in this case the anomalous and the MSSM sec- 
tors are decoupled, we can choose the model parameters 
to keep fixed the mass gap between the LSP and the NLSP. 
So we have studied the behaviour of the model in function 
of the mass of the LSP for a defined mass gap. 
We have studied separately the two cases in which the 
NLSP is mostly a bino and in which the NLSP is mostly a 
wino. In the first case the coannihilation involves only the 
LSP and the NLSP while in the latter case it involves also 
the lightest chargino, that in the MSSM is almost degener- 
ate in mass with the wino. So they are different situations. 
The results of the relic density calculations are showed in 
figures [1] and [5] for the bino NLSP case and in figures |3] 
and [3] for the wino LSP case. 
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Figure 1: Relic density for a bino-higgsino type NLSP and mass gap 
1 % 
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Figure 2: Relic density for a bino-higgsino type NLSP and mass gap 
5% 
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Figure 3: Relic density for a wino type NLSP and mass gap 5% 
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Figure 4: Relic density for a wino type NLSP and mass gap 10 % 
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156 In those figures we have shown the most significant209 

157 results, but we have found model that satisfy WMAP data2io 

158 up to 10% mass gap for the bino NLSP case and up to 20%2ii 

159 for the wino NLSP case. 212 

213 

160 4-^- Mixing case m 

161 No mixing implies that Qh^, 7^ 0, so we can't have con-^'^ 

162 trol on the mass gap between the LSP and NLSP, because^"^ 

163 the mass matrix is no more decoupled in two blocks. So we 

164 choose to let all parameters unconstrained and therefore 

165 to collect data in the mass gap ranges 5%, 5 10%. 

166 In each case we have started our study scanning the pa-219 

167 rameter space in search of the region permitted by the22o 

168 experimental and theoretical constraints, i.e. the region 

169 in which we could satisfy the WMAP data with a certain^" 

170 choice of the model parameters. After that we have nu-^^^ 

171 merically explored this regions to find sets of parameters 

172 that satisfy the WMAP data for the relic density. We have 

173 found many suitable combinations for both types of NLSP. 

174 So we have chosen some of these successful models and we225 

175 have computed the relic density keeping constant all but226 

176 two parameters and plotting the results. 227 

177 We have found regions of the parameter space in which228 

178 the WMAP data are satisfied for mass gap over 20%, but229 

179 in the following we will only show results for the regions 

180 5% and 5 10%, because they are more significant. 230 

181 For simplicity we will refer to these regions as "5% region23j 

182 "and "10% region" respectively. 232 



4-2.1. General results 235 
In this section we want to list some results that are236 
valid for both types of NLSP. 237 
First of all, given the constraints on the neutral mixing de-238 
scribed in [Sl, we have obtained that — 1 ^ Qh^ ^ 1- This239 
constraint was not appearing in our preliminary analysis24o 
in ^ . This implies that also in our general case the mixing24i 
between our anomalous LSP and the NLSP is small. 242 
We have checked that there are suitable parameter space243 
regions in which the WMAP data are satisfied and we have244 
found that this is true for all possible composition of our 
anomalous LSP. 

We have also checked that in each region we have studied^*' 
there are no divergences or unstable behaviours in our nu- 

° 246 

merical results. 

247 

We have verified that the relic density is strongly depen-^^^ 
dent on the LSP and NLSP masses and composition while 

^ 249 

it is much less dependent on the other variables. Anyway 

r 250 

it can be shown that there are cases in which the parame- 

^ 251 

ters not related to the LSP or NLSP can play an important^^^ 
role. We will show an example of this case in a forthcoming^^^ 
subsection. 



We chose two sample models which satisfy the WMAP 
data [18| with mass gaps 5% and 10%. We study these 
models to show the dependence of the relic density from 
the LSP composition and from the mass gap. To obtain 
these result, we have performed a numerical simulation 
in which we vary only the stiickelino and the primeino 
soft masses. The results are showed in figure 1, with the 
conventions: 

• Inside the continous lines we have the region in which 

{nh'^)wAiAP ^ ^h"^ 

• Inside the thick lines we have the region in which 
{Vtli^ ~ 'i(T)wMAP < < i^lh^ + 3o-) 



JWMAP 



206 4--2.S. Bino-higgsino NLSP 

207 If the NLSP is mostly a bino-higgsino we have a two, 

208 particle coannihilation. 



• Inside the dashed lines we have the region in which 
{Vtli^ - b(T)wMAP < ^h'^ < iS^h"^ + ^cf)wmap 

• Inside the dotted lines we have the region in which 
{nii^ ~ lQa)wMAP < ^h"^ < i^h'^ + 10(t)wmap 

Going from the region with a 5% mass gap to that with 
a 10% mass gap there is a large portion of the parameter 
space in which the WMAP data cannot be satisfied, while 
the regions showed in the second and third plot in fig. 1 
are similar and thus are mass-gap independent. 

4.2.3. Wino-higgsino NLSP 

If the NLSP is mostly a wino-higgsino we have a three 
particle coannihilation, because the lightest wino is almost 
degenerate in mass with the lightest chargino, so they both 
contribute to the coannihilations. 

In this case we perform the same numerical calculation 
illustrated in the previous subsection. We have extensively 
studied a sample model with mass gap 10%, showing an 
example of funnel region, a resonance that occurs when 
2 Mlsp ^ Maq. In our sample model Mlsp ^ 330 GeV, 
while Mao ~ 630 GeV and this leads to the relic density 
plot showed in figure 2, with the same conventions used in 
figure 1. So we can state that also in the case of anomalous 
LSP we can have a behaviour similar to that of the MSSM, 
given that the LSP coannihilates with a MSSM NLSP. 

5. Conclusion 

We have modified the DarkSUSY package in the rou- 
tines which calculate the cross section of a given supersym- 
metric particle (contained in the folder '-^/src/an) adding 
all the new interactions introduced by our anomalous ex- 
tension of the MSSM. We have also written new subrou- 
tines to calculate amplitudes that differ from those already 
contained in DarkSUSY. We have also modified the rou- 
tines that generate the supersymmetric model from the 
inputs, adding the parameters necessary to generate the 
MiAUSSM ^] and changing the routines that define the 
model (contained in ^/src/su) accordingly. Finally we 
have written a main program that lets the user choose if 
he wants to perform the relic density calculation in the 
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Figure 5: Plot of the relic density of the LSP vs the stiickelino and289 
primeino masses. The first plot shows the case of mass gap 5 %-29o 
The second plot is a zoom of the first in the region between 20-140 
GeV for stiickelino mass and 500-850 GeV for primeino mass. The 
third plot shows the case of mass gap 10 % in the same region of the^'^ 
second plot. All the masses are expressed in GeV 
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297 



MSSM or in the MiAUSSM. The code of our version of 
the package is available contacting andrea.mammarella 
@roma2. infn.it. 

These modifications have permitted to extensively numer- 
ically explore the parameters space for an anomalous ex-^'" 
tension of the MSSM without restriction on the neutraljoo 
mixing or on the free parameters of the model. We havesoi 
verified that our model does not lead to any divergence or 
instability. 
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Figure 6: Example of a funnel region for a wino-like NLSP, whose 
mass gap with the LSP is 10% 



ing between the anomalous and the MSSM sectors and the 
case in which there isn't such mixing. 
In the no-mixing case we are able to keep fixed the mass 
gap between the LSP and the NLSP. We have found that 
if the NLSP is mostly a bino we can satisfy the WMAP 
up to 10% mass gap; if the NLSP is mostly a wino we can 
satisfy the WMAP up to 20% mass gap. 
In the mixing case we have we have obtained a constraint 
on ■ We have also found sizable regions in which we 
can satisfy the WMAP data for mass gaps which go from 
5% to beyond 20%. We have studied some specific sets of 
parameters for the 5% and 10% mass gap regions, showing 
that relatively small changes in the mass gap can produce 
very important changes in the area of the regions which 
satisfy the experimental constraints. 

We have also showed, as an example that we still have the 
MSSM physics, the presence of a funnel region, analogous 
to that in the MSSM, in some region of the parameter 
space. 

So we can say that a model with an anomalous LSP can 
satisfy all the current experimental constraints, can show 
a phenomenology similar to that expected from a MSSM 
LSP and can be viable to explain the DM abundance with- 
out any arbitrary constraints on its parameters. 
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268 We have studied separately the case in which there is mix- 
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